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ABSTRACT 



Open systems acquire time-dependent coupling constants through interaction with 
an external field or environment. We generalize the Lewis- Riesenfeld invariant theorem 
to open system of quantum fields after second quantization. The generalized invariants 
and thereby the quantum evolution are found explicitly for time- dependent quadratic 
fermionic systems. The pair production of fermions is computed and other physical 
implications are discussed. 
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The dynamical evolution of the Universe in time has to be examined carefully to understand 
the phase transitions that are crucial to the cosmological scenarios. The phase transition may 
be affecting topologically unstable structures in the Universe and these may have important 
consequences for the inflationary Universe. The observational aspects like baryon asymmetry 
may require the presence of time-dependent gauge fields. The time evolution of the Universe 
would require that both scalar and fermionic fields are time dependent through their interactions 
with the gravity. Furthermore one may need to calculate the pair production of fermions due to 
the oscillating inflatons. This may provide a possible mechanism for preheating mechanism for 
fermionic systems. 

All these time-dependent processes in the Universe would be happening when the system is not 
in equilibrium and there is no fixed temperature. We have a system that is time-dependent and in 
a non-equilibrium state. Similar time-dependent and rapidly varying processes in the laboratory, 
that simulate the conditions at the early evolutionary stage of the Universe, relate to the heavy 
ion reactions at the relativistic heavy ion facility at Brookhaven. It is intended that in colliding 
heavy ion (U on U) reactions there will be a deconfining phase transition from nucleonic to a 
quark-gluon system that exists for a very brief period of time before making a phase transition 
back to the baryonic matter. The time evolution of the formation and dissolution of the quark- 
gluon plasma is a highly non-equilibrium process that requires a careful consideration of the time 
dependence of the coupling of both the fermionic and bosonic systems. 

The coupling constants or parameters of a system depend on time implicitly through the 
interaction when it interacts with an environment. Any system whose coupling constants or 
parameters depend on time explicitly can be regarded as a subsystem of a larger closed system 
and belongs to an open system. To describe properly the quantum states of such open system one 
may find the evolution of the system plus environment and integrate out the degrees of freedom 
corresponding to the environment. Or one may rely on field theoretical methods such as closed 
time path integral (for review and references, see [^) or thermo field dynamics (for review and 
references, see 0). 

Even though the thermo field dynamics is a real time operator formalism it has not been 
applied to a real problem to investigate the time evolution. Hence we choose instead to use the 
functional Schrodinger picture approach [3-5]. It is conceivable that this may provide some help 
in using the thermo field dynamics. Even the quantum evolution of a time- dependent system has 
also been an outstanding problem in quantum mechanics. A frequently used method is to find 
trial wave functions (a variational approach) with undetermined time-dependent coefficients and 
to fix these coefficients by solving the Schrodinger equation. Lewis and Riesenfeld introduced a 
systematic method to find the operators that satisfy the quantum Liouville-von Neumann equation 
and to solve for the exact quantum states as their eigenstates 0. This method has been applied 
successfully to some of the time-dependent systems, in particular, the harmonic oscillator systems. 

The Lewis-Riesenfeld invariant theorem has been applied successfully to study the quantum 
evolution of scalar fields in an expanding Friedmann- Robertson- Walker universe [7-12]. A variety 
of results were obtained including the vacuum state and the Wightman functions with respect to 
the vacuum, thermal equilibrium and coherent states. However, as suggested above the overall 



evolution of the universe would require not only the time dependence of the scalar field but also 
that of a fermion field. Since the fermion field is a four component spinor, it appears natural to 
use the Lewis-Riesenfeld invariant theorem to the fermion field taking proper account of the fact 
that the spinor components lead to a set of coupled equations. The usual procedure of writing the 
spinor field in terms of creation and annihilation operators is adopted. The results of quantum 
evolution of the fermion field are presented. We give a brief overview of the application to the 
scalar field and then a detailed account of the method as applied to the fermion system is given. 

The open system of the quantum field is described by the functional Schrodinger equation 

[3-5] 

z^|vl>(0,t)) = i^(0,-zA t)|vl>(0,t)), (1) 

where (j) represents a scalar or fermion field. Our stratagem is to express the Hamiltonian in the 
second quantization as a sum of infinite number of finite dimensional quantum systems, which 
are decoupled or coupled depending on the type of interactions. And then we apply the Lewis- 
Riesenfeld theorem to each finite dimensional system and find the generalized invariants and 
thereby the evolution of the second quantized Hamiltonian. 

According to the Lewis-Riesenfeld invariant theorem, the hermitian invariant operator 0{t), 
satisfying the Liouville-von Neumann equation for a time-dependent Hamiltonian H{t) 

^|^W = 4^W + [0(t),i/(t)] = O, (2) 

has time-dependent eigenstates and time-independent eigenvalues 

d{t)\Xn,t) = X„\Xn,t). (3) 

And the exact quantum evolution of the system 

z^^Mt)) = mmt)), (4) 

is given by a superposition of states |A„,t), up to time- dependent phase factors. 

The theorem was originally proved by Lewis and Riesenfeld and applied to a time-dependent 
harmonic oscillator |^. We now give a simple proof of the theorem. We assume that 0{t) be a 
hermitian operator and act on the Hilbert space for the system H{t). When the states are not 
normalizable, the Hilbert space should be extended to a rigged one [|l3l. The operator may have 
a possible degeneracy 

dit)\xit) = x^\x:,t), (5) 

where the index a runs over the degenerate eigenstates. One can apply eq. (||) to the eigenstate 
|A5^,t) and differentiate eq. (H) to obtain the following identity 

where / denotes the identity operator. By acting (A°, t| on eq. (|) and using the fact that 0{t) is 
the hermitian operator and (A^, t\ belongs to the kernel of the hermitian operator 0{t) — Xnl, one 



finds that ^^ = 0, that is, the eigenvalue is a constant. Hence the first term of eq. (^) vanishes 
and eq. (||) reduces to 

{^^^- m)\K,t) = T.^nit)\>^lt). (7) 

Assuming the orthonormal basis (A5^,t|A^,t) = Smn^ab, f'^i.t) is determined by acting (A[^,t| on 
the both sides: 

v':it) = {xit\{i^^- m)\xit). (8) 

Since 0{t) — Xnl is a hermitian operator and ff^it) is an element of the corresponding hermitian 
matrix, there exist a set of orthonormal eigenvectors such that 

^'::{t)Cait) = ^:itmt). (9) 

For the new basis for the subspace spanned by eigenvectors |A^, t) for a given A„: 

K,t) = J2Cait)\xit), (10) 

a 

one then has 

{i^^-H{t))\x:,t) = ^:{t)K,t), (11) 

from which follows the equation 

V^it) = {K.t\{^§-^- m)\Xlt). (12) 

The quantum state given by 

|vI/„,t) = e^/<W|A^,t), (13) 

is an exact solution of eq. (^. The most general state is a superposition of states like the ones 
given by eq. (|TBp. 

It should be pointed out that the Lewis-Riesenfeld theorem can also be applied to infinite- 
dimensional quantum systems such as the quantum fields when they are expressed as a sum 
of finite dimensional systems. The second quantized Hamiltonian of quantum fields has such 
a form. Another important observation is that the theorem can also be applied to bosonic or 
fermionic systems when the commutators or anticommutators are appropriately used for boson 
and fermion operators, respectively. We shall consider a scalar field with a time-dependent mass 
and a time-dependent fermion system with the most general quadratic form as open systems. 



A massive scalar field with a time-dependent mass described by the Lagrangian density 

C=^-g^^d,cpd^<l>-'^<P' (14) 

has the second quantized Hamiltonian in the standard basis of the Minkowski spacetime 

m = E^i^'HOK^^ + l)+ ^o.{t){ai' + al) ^EHait). (15) 



Here, a denotes collectively the Fourier sine/cosine modes, and 

fii"Ht) = f + ^, nut) = ^-\, (16) 

where 

ul = m^it) + k^ (17) 

The creation and annihilation operators for each mode satisfy the usual commutation relation 

[ac„al] = 5ap- (18) 

The Lewis-Riesenfeld theorem applied to Ha{t) gives rise to a pair of invariants [7-12] 

Mt) = fi'\t)a^ + f^o^\t)al Aiit) = h.c, (19) 

where 

fL-\t) = ^i-^^:{t) + vut)), e\t) = ^H^iit) - Km (20) 

(^«(t)+a;2(t)(^„(t) = 0. (21) 

It should be remarked that V5k(^) satisfies the classical equation for the k-mode. The equal time 
commutator 

K(t),4(t)] = 5,^, (22) 

is satisfied for all times for a complex solution due to the wronskian of eq. ( pi]) : 

t{Mi)¥^*ait) - KMt)) = ^- (23) 

The number operators and their eigenstates 

Na{t) = Ai{t)A^{t), N^K,t)=n^\n^,t), (24) 

define the Fock space for each mode. The exact quantum states for each mode are linear super- 
position of these number states as given by eq. (|13D, and the wave function for the scalar field is 
the product of these quantum states. 

We now generalize the quadratic form of eq. ([T5|) to a fermionic system by defining the 
Hamiltonian 

a 

+Qi-+\t)bJi + Ql+-\t)bid^. (25) 

Here, a denotes a collective notation for the momentum k and spinor index r, and X^a = 
/(i^kX)r=i2 5 S'^d the particle and antiparticle creation and annihilation operators satisfy the 
anticommutators, respectively, 

{ba, bi} = Sad, {da, dU = S^p, (26) 



and all the other anticoinmutators vanish. So the Hamiltonian given by eq. ( pSj ) describes a 
system bilinear in the fermion fields. The unitarity of quantum evolution requires the hermiticity 
of the Hamiltonian which leads to the conditions 

It is found that an invariant of the form 

Oa{t) = fL-\t% + f^^\t)hl + gi-\t)L + 9^^\t)dl (28) 

satisfies eq. (|) provided that 

4^<^"^(^) + ^i''H^)/^nt) - fii-+Ht)^(-)(t) - fiL-)(t)^W(t) = 0, 
^|/^+Ht) - fiL^Ht)/W(t) - ^^^^\t)gV{t) - nL^-\t)gi^\t) = 0, 
^^/a~\t) + ni^\t)gi-\t) + n^^-\t)f!f\t) + nl-\t)fi^\t) = 0, 
^l^i+H^) - ^lF\t)gM(t) + ni^^\t)ft\t) + nl-^\t)f!^-^\t) = 0. (29) 

By introducing two column vectors 

rnt)- ^ fftKt) + f!.-'Kt)\ 1 fgi-Kt)+9k^Kt)\ .on^ 

^"^'^ - 71 \fi-\t) - /w(t) j ' ^-^'^ - 71 Ui-)(o - 9i^m ) ' ^3°) 

we can write eq. ( pQ]) as two vector equations 

^^^a(t) + nl''\t)a,Ua{t) + MuSt)Va{t) = 0, (31) 

z^y„(t) + fi(f )(t)aiK(t) + MvAt)Ua{t) = 0. (32) 

Here, Mu^{t) and My^(t) denote the mixing matrices between Ua{t) and Va{t): 

MuAt) = AL°)(0/ + A«(0^i + A(^)(t)a2 + Ai3)(t)a3, 
Mv.(t) = -A(°)(t)/ + A«(t)ai - A(^)(t)a2 - A(f)(t)a3, (33) 

where (Tj are Pauli spin matrices and 

AL°Ht) = -l{ni-^\t)-ni~^>it)), A«(t) = -^K-+)(t) + !](-+> (t)), 

Ai^Ht) = l{^L-\t)+^L-'*it)), A(f)(t) = i(fi(++)(t) -fiL++>(t)). (34) 

By eliminating either Ua{t) or Va{t) in eq. ( p2D we obtain the second order differential equations. 
So we get two types of annihilation operators 

Bait) = fi^-\t)b^ + /f +)(t)6i + 9L''-\t)d^ + 9i'''-\t)dl 

Da{t) = ff-\t)K + /f +)(t)^l + 9?-^it)L + 9i''''\t)dl (35) 



where /^^ \t) and g^J^ \t) are chosen the (adiabatic) positive frequency solutions. 

One can make Ba{t),Da{t) and their hermitian conjugates the annihilation and creation op- 
erators that satisfy the anticommutators at each time 

{B^{t),Blit)} = 5a(,, {D^it),Dlit)} = 5a(3. (36) 

These lead to the consistency condition for the coefficient functions 

/f-^*W/f-nt) +/f +)*(t)/f +nt) +^i^-)*(t)^i^-)(^) +^i^+)*(t)^i^+)W = 1, 
e-^%t)fr\t) + fr^*it)e^\t) + 9r>{t)g^f-\t) + 9i''^>{t)gi^^\t) = 1. (37) 

The consistency of anticommutators for all times is guaranteed by 

|-(f/(^)t(t)^(B)(^) + \/(^)t(t)y(^)(t)) = 0, 
^^{lJ^^)^{t)U^^\t) + y(^)t(t)y(^)(t)) = 0, (38) 

where M^^^t) = Mv„(t) is used. The other anticommutators 

{B^it),Bpit)} = {BUt),Bl{t)} = 0, 

{D„(t), D^(t)} = {Diit),Dlit)} = 0, (39) 

are consistently satisfied for all times due to 

|(^(^)t(,)^3^(B)(,) + y(^)t(t)a3\/(^)(t)) = 0, 

|([/(^)t(i)^3^(D)(^) + V(^)\t)asV^^\t)) = 0, (40) 

where (T3M^^(t)o'3 = —My^^t) is used. Still other anticommutators 

{B,{t),Dp{t)} = {BUt),mt)} = 0, 

{B^{t),Dl{t)} = {BUt),bfs{t)} = 0, (41) 

are also satisfied for all times due to 

|([/(^)^(t)a3t/(^)(t) + V^^^^it)a,V^^\t)) = 0, 

^^{u^^)\t)U^^)(t) + y(^»(t)\/i^)(t)) = 0. (42) 

Since the Ba{t),Bl^(t) and D ait) , D^^it) share the same properties of the annihilation and 
creation operators for time-independent system, they will be used as those for time-dependent 
system. So the number operators are defined by 

N^^\t) = Bi{t)B.{t\ N^^\t) = DUt)Da{t). (43) 



The eigenstates for each mode are the number states of NJ^^\t) or Nj^^\t). For instance, the 
zero-particle state for NJ^^^t) is given by 

|0(f\t) = i?(t)|0L^0W) + cL?(t)|0(^l('^)) +ei?(t)|l(^0('^)) +ei?(t)|l^ (44) 

where \n'^^\n^'^^) are the number states with respect to bj^ba and dl^da, respectively, and 

f?-\t)e^\t) - e-\t)e^\t) = 0. (45) 

One-particle state is constructed by applying the operator Bj^lt) to the zero-particle state (eq. 
(0)). The number states of the antiparticle can be similarly constructed by using the number 
operator N^^\t). These states in general depend on time, but their eigenvalues (occupation 
numbers) do not depend on time as expected according to the Lewis-Riesenfeld theorem. The 
vacuum state is the one that is annihilated by all the annihilation operators Ba{t) and Da{t): 

S,(t)|0,0,t) = D,(t)|0,0,t) = 0. (46) 

When the interaction is turned on for a finite period of time or the system evolves from an 
asymptotic region where all the Via are constants and eq. (^61) coincides with the Minkowski 
vacuum, the pair production of fermions from the initial vacuum is given by 

(o,o| Y.i^'J'Kt) + N!^\t))%o) = E(/f ^^*W/f ^^ + ^i''^^*W^f ^0- (47) 



Let us consider some particular examples of using the Hamiltonian given by eq. (^). 

Dirac Oscillator. It is the case where i^^ ¥" ^^d i^^ — 0- Then the first order equations (eq. 
(P9|)) can be easily integrated as Ba{t) = e^J " b^ and Da(t) = e*J " da and their hermitian 
conjugates. So the generalized invariants become B]^{t)Ba{t) = b\ba and D\^{t)Da{t) = d\da- 
Their number states are not only the eigenstates of the generalized invariants but also those that 
diagonalize the Hamiltonian. 

Fermion Scalar Field Coupling. The case where f^^""^-* = ^^^ ■* — 0? ^iF^ = ^0 + ^2! ^^^ 
Qi—) = -Q(++) = ^^:^ leads to the Hamiltonian 

mg(j){t) 



H{t)=Jd'k {ko + ^^) Y. {bli^M^) - l{k)dl{k)} 



r=l,2 
t) , 

2ko ^^12 



mg(j){t) 



Y, {bl{k)dl{k) + dr{k)br{k)} 



(4J 



This describes a fermion coupled to a uniform scalar field through the Yukawa coupling 

C = i(^^,^M(a^^) _ (a^^)i7'^V') - m^tfj - \(t>m^, (49) 



where A is a coupling constant. The coefficient functions are found to satisfy 



+{o™^ + n(-)^ + in>->l-{%^)}{ft^ + gM) = 0. 









The positive solution fj^~\t) and the negative one 5'*-~^''(t) with the condition (^) give the correct 
pair production rate. Interestingly enough the free fermion in the Friedniann-Robertson- Walker 
universe has the same form of Hamiltonian as eq. ( ^H]) |T^ . 



A few comments are in order. The generalized invariants in eq. ( p5D for the Hamiltonian 
given by eq. ( P5| ) are more general than those for eq. (^5D, since Ba{t) involves only b^ and 
dl^ and Da(i) involve only b^, and da- When one has cross terms like bad]^ and bl^da, then eq. 
(^) is the unique choice compatible with the algebra. Furthermore, eq. (47), when restricted to 



eq. ( ^81) gives the same result for the pair production as other quantum field theories. One may 
also extend the analysis to the fermion interacting with an electromagnetic field by introducing a 
Hamiltonian that is similar in form to eq. (|25| ) but involves Dirac matrices 7^ in a coupling with 
the time dependent coefficient fiait)- 

One of the important issues not treated in this paper is the non-equilibrium aspect of open 
system of quantum fields. This can be done by introducing a density operator of the form 
Pa{t) = e~^*^^" (t)+i^a (t))^ which satisfies the Liouville-von Neumann equation by its construction, 
which will be addressed in a future publication. One of the applications of the result of this paper is 
the evolution of fermions in an expanding Universe, where fermions gain time-dependence through 
the interaction with the gravity. The other is to calculate the pair production rate of fermions 
due to the oscillating infiaton and to provide a preheating mechanism for the fermionic system. 
Finally the time evolution of the quark-gluon plasma in a relativistic heavy ion collision would 
simulate the conditions present at the early moments of the Universe and would thus provide us 
with an experimental view of the evolution of the Universe. 



After the submission of this paper we were informed of Ref. [jT5[, where linear invariants 



have also been introduced for fermionic systems. However, we have explicitly shown the auxiliary 
equations for invariants, proved the consistency of the method and constructed the Hilbert space 
of particle and antiparticle states for general fermionic systems. 
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